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Предлаг ается численно-аналитический метод решения несмешанной осесимметричной задачи термоупруго- 
сти для непрерывно-неоднородного по Г; лубине полупространства при заданных на его поверхности источни- 
ках тепла и силовом воздействии. Для решения г! раничной задачи используется аппарат интегральных пре- 
образований Ханкеля. Приводится анализ влияния изменяющихся по Г. лубине Термомеханических парамет- 
ров в приповерхностном слое на распределение температурного поля, теплового потока и смещение по- 
верхности. 

Ключевые слова: термоупругость, непрерывно-неоднородное полупространство, численно-аналитический 
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1. Постановка граничной квазистатической задачи термоупругости для непрерывно- 
неоднородного по глубине полупространства при заданных на его поверхности усили- 
ях и источниках. Рассмотрим полупространство, упругие характеристики которого непрерывно 
меняются с глубиной в пределах прилегающего к поверхности слоя толщиной #Н, а затем стабили- 
зируются и остаются постоянными (рис. 1). С полупространством свяжем цилиндрическую систе- 


му координат (7,ф,2). Обозначим через и,у, » смещения вдоль осей уф, 2, ные 
радиальное, угловое, нормальное и тангенциальные напряжения соответственно. Кроме того, 


введем обозначения для температуры - 7, коэффициента теплопроводности - ^.,(2), коэффици- 
ента линейного расширения -— «а, (2), коэффициента теплоемкости -— с,(2). В ненапряженном со- 
стоянии температуру полупространства примем равной Т, и будем анализировать разность тем- 
ператур 9=Т-Т.. 








Р(г) 0’(г) 
=5_^ 
т 
а 
р Е (=) № (=) 
А, (=) @ (=) 
|. #2 





Рис.1. Схема нагружения 
непрерывно-неоднородного полупространства 
Для определения полей перемещений, деформаций, напряжений и температуры в рамках 
осесимметричной квазистатической несвязанной задачи термоупругости мы имеем следующие 
уравнения: 





* Работа выполнена при финансовой поддержке грантов РФФИ (09-08-011410а, 10-08-01296-а, 10-08-90025-Ве!-а), 
ГК № 02.740.11.0413, ГК № 02.740.11.5193, ГК № Р1107, АВЦП 2.1.2/10063. 
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1. Уравнения равновесия неоднородного по глубине полупространства при отсутствии 
массовых сил, записанные в цилиндрической системе координат: 








д От т 

1, 2", +" 0 ео (1.1) 
д Оба. 
э=(”т,. |+ 02 =. 

2. в Е (1.2) 


2. Уравнение неоднородной по глубине теплопроводности для о поля: 
2 2 2 д), 
ов о + т (2) 00 





а = 92 5+ ст (2) 57 =0. (1.3) 
3. Уравнения связи термоупругих напряжений с деформациями (Дюамеля-Неймана) [1]: 
0, =2М(2)Е, +Л(2)8-3К (2). (2)0, 
<,= 2М(2)= ты Л(2)=-3К (2). (2)0, 
02=2М(2)=, +Л(2)5-3К(2)%.(2)9, 
т 2М(2) Е, = МЕ, (1.4) 
= 2М(2)= (=› 
ЗК(2)=3/Л(2)+2М (2), 
9=7Т-7х, &=8, +8, +8.. 


Выражения для компонент деформаций имеют вид: 











ди и ди’ 
И 
г № фт 2 д 
. _1 ди ди Е этой - _ Иду у (1.5) 
2 2| 02 ди $ 2 02 5 2 б-в 
ай 
бу г 02. 


Далее предполагается, что как коэффициенты Ламе М(2) и Л(2), так и коэффициенты 
теплопроводности Ат(2), линейного расширения 0, (2), а также теплоемкости С,(2) являются 
непрерывными функциями координаты 2, такими, что 


1. М(2)=М(-Н)=М*у, —®<2<-Н, 
А(2)=АСН)=Л», ^,(2)=№(-Н)=Ат; 
0, (2) = ©, (-Н)=ох. 
2. М(2)=МС(2), А(2)=ЛС(2), -Н<250, (1.6) 
^.(2)=№:(2), @.(2)=95(2), с,(2)=ст(2). 
3. М“(-Н)=М*, Лб(-Н)=А5, 
№(-Н)=Ал. о (-Н)=@х, сг(2)=ст(2), 
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где Н -— толщина неоднородного слоя, сцепленного с подстилающим полупространством, то есть 
глубина, с которой мы полагаем термомеханические характеристики полупространства по- 
стоянными. Индекс $ соответствует подстилающему однородному полупространству, а С- 


* * * * * 
неоднородному слою; Л, Л, М, М, А.., Л... а.., @., с,., с; - произвольные кон- 
* * 


станты. 

Наряду с парой коэффициентов Ламе для описания упругого поведения твердого изотроп- 
ного тела используются: модуль сдвига С и коэффициент Пуассона У или модуль Юнга Е и ко- 
эффициент Пуассона У. Коэффициенты Ламе Л и М (иногда обозначаемый С и называемый 
модулем сдвига) связаны с модулем Юнга Е, коэффициентом Пуассона У и изотермическим моду- 
лем объемного расширения Ксоотношениями: 


д, о Еу _ М(2М+3^) 
Му А+) "= МА ти 
у. ЗК) ЕЗА М) 


2(М+лЛ)” 

Запишем на основании (1.4), (1.5) представления для напряжений через смещения и под- 
ставим их в уравнения (1.1) - (1.2), получаем систему дифференциальных уравнений в частных 
производных второго порядка относительно перемещений и температуры (уравнения Ламе). Эти 
уравнения могут быть представлены в форме 


2 ди» ди 
М(=) у мени) МС |2" 5; и 
00 
= (2) >, 
меч (ме)+л (2) -+2М “(2 АЕ (1.8) 
=^’(2)0+^(2)0', 


5 У р ду т 
М(2) У +М о 


Г 


^.,(2)У20-+^, (2)0'= 





Здесь используются следующие обозначения: 


Мк2= АМ(2) ии 4А(2) 72-19|,9 |, 6? 
42 42 





д’ 








Г ОГ д22 ’ 
ИО, а - 4, =, и -®. 
д г 02 4 9 е 


Мы полагаем, что поверхность полупространства нагревается в пределах круга радиусом а 
непрерывно действующим тепловым потоком. Вне области нагрева поверхность полупространства 
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теплоизолирована. Предположим также, что внутри круга приложена произвольная нормальная и 
касательная нагрузка. 


0'(",0)=-аВ(г), г<а, 


0'(",0)=0, г>а, 

с,(г,0)=-р(г), 0<г<а, 2=0, 

о.(7,0)=0, аи. 2-0, (1.9) 
с. (г,0)=-—а(г), 0<г<а, 2=0, 

с: (7,0) =0, а<г<®, 2=0, 

т. (7,0) =Кг), 0<и<а, 2=0, 

т., (7,0) =0, ага, ==0, 


Здесь а — коэффициент теплопроводности. 

На границе сцепления неоднородного слоя с однородным полупространством при 2=-Н 
в силу непрерывности должны выполняться условия сопряжения по смещениям, напряжениям, 
температуре и тепловому потоку. 


о (’-Н) — а, ("-Н), Е ("-Н) = Е. (“-Н), 


иб("-Н)=и° ("-Н), С (и-Н) = ("-Н), (1.10) 
09° ("-Н)=0° (,.-Н),(0°)'(=-Н) = (05) (*,-Н). 


На бесконечности при (’,—2)— ® смещения, деформации и напряжения исчезают. Зна- 
чения разности температур и теплового потока при этом также стремятся к нулю. 


Пит | №,\,5.,8 ,5 ‚5 ‚© о бот ,0,0'’=0 
7 — 00 э 2 и? ф’ = 12? ры ф’ 7? тд? 2 ь (1.11) 
. \_ 

[име ее в 6,0,б от 0,0 )- 0. 


Таким образом, мы сформулировали первую граничную квазистатическую осесимметрич- 
ную термоупругую задачу для неоднородного полупространства: найти смещения, деформации, 
напряжения и распределение температуры внутри полупространства, удовлетворяющие системе 
дифференциальных уравнений (1.8) при заданных граничных условиях и распределении воздей- 
ствий на границе (1.9) - (1.11). 

2. Построение фундаментального решения квазистатической осесимметричной термо- 
упругой задачи для неоднородного по глубине полупространства. Будем разыскивать 


решение для смещений ци, уи и/и температуры © в виде интегралов Ханкеля [2]: 


и (г.г) =, 2) (иди, и.) = ТУ) рум, 
0 0 


[е.0) [е.0) 
и{г,2)= [И (1,2) (ууу, 0(т,2)= 1 Т(р2) 7 укуиау. 
0 0 
Подставим (2.1) в систему дифференциальных уравнений в частных производных (1.8) и, 
приравняв к нулю подынтегральные выражения, получим систему обыкновенных дифференци- 
альных уравнений второго порядка 


(2.1) 
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МИ" +у(М+ МИ '- у? 22М+^А)0 +М'О'+ УМ =АГ, 

ОМ+лА)И”" — (М + МИ '- УМ +2М'+ Л), - УЛ = 
=АТ +АТ’,, 

МУ" + М'У' -у*МУ =0, 

де (2 (ТР) - у е ТС +“ = 0. 


(2.2) 





(здесь ' указывает на дифференцирование по 2). 
В системе дифференциальных уравнений (2.2) третье уравнение не связано с остальными, 


и его решение здесь рассматривать не будем. 
Граничные условия (1.9) примут следующий вид: 


(2м(0)+^(0)} и.) -УА (00,0) ()Т(у,0)=-Р(, 


м(0) (у (у,0)+0 (4.0) =О(у, 


Т'(у.0) ==@В(У), (2.3) 


Рд= ТРУ фурар, 96) = ау ФОрар. 
0 0 


[9.9] 
В(у)= | ВФУУ (рурар. 
0 
Используя векторное представление для трансформант 
Е 
хо =(х,х, ХХ), 
х=ИЦ, х,=И’, х,=Й’, м. =Й”, (2.4) 
Г 
запишем систему (2.2) с исключенным третьим уравнением в матричном виде, при этом явно вы- 
делим части, соответствующие покрытию и подложке: 


ахС 




















^^ =АСхе, -Н<2<50. (2.5) 
Г 
0 1 0 0 0 0 
22М+А мМ' ты М+л ук 0 
М М М М М 
0 0 0 1 0 0 
АС = Л М+л ›‚ М 2М'+Л' Г к |: 
МЛ МА ' 2М+л 2М+^ 2М+л 2М+Л 
0 0 0 0 0 1 
с’ 
0 0 0 0 у? = 
у 
5 
А, -®<2<-Н. (2.6) 
И 
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0 1 0 0 0 0 
22М+А 0 Их У 0 
М М М 
д 0 0 0 1 0 0 
0 Е у? М 0 К 
2М+А 2М+А 2М+лА 
0 0 0 0 0 1 
0 0 0 0 у’ 0 


Граничные условия при этом имеют вид: 
(2м(0)+А (0) (у,0)-уЛ(0)хг (у,0)-К(0)х; =-Р(у), 
МО (ноухс 6 =0(, (27) 
х (1,0) == В(); 
хС (у,0) = х(у,-Н), хС (у,0) = х(у,-Н), 
1 1 3 3 
хб (у,0) = х° (у,-Н), хС (у,0) = х°(у,-Н), 
2 2 4 4 


хС (у,0)=х5 (у,-Н), хС (у,0) =х°(у,-Н). 
5 5 6 6 


Необходимо отметить, что последнее условие имеет место в силу непрерывности измене- 
ния свойств слоя и основания (1.5). 
Общее решение системы (2.2) для однородного полупространства 


А'= М'=К"=А, =0, М>0, Л>0, К>О0, ^.>0 имеет вид: 
$ я г 
х8(у,2)= (4 +2, + (к, +к,2+ ка, с | 
х(ь2) =(а, + +24, + (к, Гук + кии), уе", 


о (,2) = | ЕК у2)4, + |-5 +к,2 + ке =" (2.8) 


х5(ь2) = (а +@-к, +124, + (к, +куг+ куда, уе", 


(уе ее", ое еГ, 


й 2 =——_ к, Е 
АМ 202М+Л) 2м 


где 4, (1=1,2,3) — произвольная функция параметра \. 


_ А+3М А р 





) 


в] 
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Решение хС(у,2) системы дифференциальных уравнений (2.5) строится методом модули- 


рующих функций. Будем искать хС(у, 2) в виде 





3 
х° (9,2) = У. А(уда (у, 2)е. (2.9) 
1=1 
Векторы а;(7,2), (1 =1,2,3) определяются из решения следующей задачи Коши: 
аа. 
: =АСа —уа, ‚изо: ре: (2.10) 
ЧР 1 1 


При начальных условиях для 2 =—Н 








а, (7,2) Ре. = (1,у,1,у,0,0) ‚ 
а, (7,2) _ „= (+12, +75у ку + 972,00); 
а, (7,2) Е (-к,+к,2+к,]2,к./у+к,2+ К), 





КЗ 
——+к,2+к,]2,к, +К,2 +к,)2,1,у) 





2=Н 


Константы 4) (: =1,2,3) определяются из условия (2.7). Таким образом, мы имеем: 


а, (у) М(у) =-Р(), 
1 


4.(у)М (у) = ОС), 


3 


= 


и [Ме 


= 


а,(у)ОКу) = -аВ(у); (2.11) 
1 


и [Ме 


1 


№) =-Л(@)уа, (1,0) + (А(0)+2М(0))а” (у,0) - (б)а> (1,0), 
Му) =М(0а, (1,0) + М(О)уа” (1,0), 
О,(у) =а° (1,0), Оу) =0,(1) =0, 
где а, (у,2), ((=1,2,3) обозначает К-ю компоненту вектора а (7, т АЕ 


&=1,2,3,4,5,6. 
а, =—№ (у)Р()- © (00) + © Фа, В(У), 


4, = ®, (у)Р(у) + ® ()0(у) + ®()аВ(у), (2.12) 
а, = ей 
О; (у) 
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м, ге О АО 9 
— к С Е С, 0,0)’ 
(М, МАЮ 1 
Е о. А =. РЕ те 0,0’ 
д, (у) = М,(у)М (у) - М (УМ У). 
Окончательно получаем следующее выражение для компонент вектора решения хС (у, =) 
при 2>2-Н: 
ж =И(у,2) = ([(у,2)Р(у) + Г.(у,2)0(у) + (у, 2)аВ(у))е" / у, 
х, =0'(у,2) = ([(у,2)Р(у)+ 1, (9,2) О(У) + 15(у,2)аВ(уе" / у, 
х; = (1,2) = (Ё(у,2)Р(у) + [(1,2)0() + В @,2)аВ(у)е" /у, 
ж =Й”'(у,2) = ([1(9,2)Р(у) + 15(у,2)О(у) + 13(у,2)аВ(у))е" / у, 
х; =Т(у,2) = [5(у,2)0В(у)е" / у, 
ж =Т'(у,2) = [8(,2)аВ(уе" у, 
[1 (9,2) = (в. (у)а>(у,2)- (ау, 2), 
[. (7,2) = (®, (у)а>(у,2) — К (у)а!(у,2))у, 


(‚2 = да, ва, 09), 
О, (у) 


№ (у) = 














(2.13) 


Введем обозначения: 


1,2) = [ Вбр, з)Ррее 1 (ку =1,21= 0, =3,4 1=0), 
0 


Г, (,2) = [У (у,2)0(ре" 1, (ук) ау, (Е =1.21=1),(& =3,41=0), 


1,2) = [14 )аВбде" 3 (указ, (Е =1,2 =, =3,4,5,6 #=0), а 
| 


Ли(г,2) = № (у,2)Р(де" 7 (к), =21=0),(К=4 1=1), 


Л.) = [Ба ,)бре" 1. (ука, =21=0)%=4 1=1), 


Ле.) = [| Ва,давбре” 1 (уе), =21= ОК =4 11) 


В соответствии с (1.4) - (1.5) мы можем теперь выписать выражения для смещений, тем- 
пературы, теплового потока и деформаций: 
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а >, 2), и(г,2) = о а (г, 2), 
и 2) 


Фи.) = ВР 16,2), 
и. 7 1, 2), 
02 ри] 02 2.15) 
" 2 8) .) 7,9-0,2) 
[@ ри Г 
Рив, = пох ле) 
[д ри 2 ри 
3 13 13 3 
== У Ль(,2)--У1(,2)+-У 1 (г,2)+У. 14 (0,2), 
1=1 Г! Е 11 
а также для напряжений: 
с, = (А(2)+2М [$ Е 1.2) + 
1=1 И 1=1 
53 з 
+^(2[8 Е, ви. 2) Е (Ры 
И 1=1 1=1 
в, = (^(2)+2М (2) У (и, =) + 
И 1=1 
3 13 3 
Ао [$ о ТИ, 10-2) = 
Ри = 1=1 
= (2) 10 (7,2), ев 


о, = (Л(2)+2М (=). 14 (7,2) + 


+2 [$ о а 1.2) ь 
1=1 И 1=1 И 1=1 

Е (а 

5. =М (2). 14(0, 2) - 2 (0,2), 


К (2) =ЗК (2)а. (2). 

3. Численный анализ решения осесимметричной граничной задачи квазистатической 
теплопроводности и термоупругости для функционально-градиентного покрытия при 
заданном на его поверхности постоянном тепловом потоке. Рассмотрим распределение 
температуры и теплового потока в покрытии, вызванное воздействием равномерного теплового 
потока с поверхности в пределах круга радиусом а, для некоторых характерных видов неодно- 
родности. 

Будем считать, что коэффициент теплопроводности /^.,(2) непрерывно изменяется в при- 
поверхностной зоне так, что выполняются соотношения (1.4). В записи граничных условий на по- 
верхности неоднородного полупространства (1.9) полагаем, что усилия отсутствуют, то есть 
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р) =а()=Кг)=0, 0<г<о, 2=0, и граничные условия на поверхности и в зоне сцепления 
неоднородного слоя с подстилающим однородным полупространством имеют вид: 
(0С)'(+,0) =-аВ(), г<а; (06)(+=,0)=0, г>а, 
0°(",-Н)=0°(х,-Н); (0°)’(",-Н) = (0) (",-Н). 
На бесконечности при (г,—2)— значения разности температур и теплового потока 


стремятся к нулю, и выполняются соотношения (1.11). 

В пункте 2 построено общее решение квазистатической осесимметричной термоупругой 
задачи для неоднородного полупространства. В развернутом виде формулы для определения 
распределения температуры и теплового потока в функционально-градиентном покрытии при за- 
данном на поверхности постоянном тепловом потоке можно записать так: 


0(",2) = а [Вора Р] 1,2)" (ру), (у")уау, 
0 


(3.1) 


(3.2) 
00(,2) Е ®° тб 2 
о В(р)рар]|› 15(у,2)е" 3, (ру) (ут ма. 
0 
Рассмотрим шесть характерных законов изменения коэффициента теплопроводности 
^.(2) в приповерхностном слое (относительно подстилающего однородного полупространства) 


(рис. 2): 
Л5(2=№.1(2), -Н<2<0, 


№: (2) = 
у 2 = соп5Е, —0<2<-Н, 
1 р | 
2)=2, )=—, 2)=2+—, Е (3.3) 
Л) №) 5 (2) Н Л (2) Е 
: а 
(2) =1-— эл / Н), т, 8 
Е(=) 
0.0 а 00 - 0.0 ВО > : 7 
0.5 в. 6 1.0 0 
1.0. ре 





Рис. 2. Схематичное изображение шести видов неоднородности 
и однородного полупространства 


На рис. 3,4 для законов (3.3) представлено изменение вдоль оси симметрии г =0 устано- 
вившегося температурного поля и теплового потока для однородного полупространства. 

При изменении коэффициента теплопроводности в 2 раза по сравнению с однородным 
полупространством температура на поверхности изменяется примерно в 1,5 раза. Видно, что из- 
менение коэффициента теплопроводности влияет только на скорость изменения температуры в 
теле, однако существенно сказывается на распределении теплового потока в приповерхностном 
слое. 
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Рис. 3. Изменение температуры для шести видов Рис. 4. Изменение теплового потока для шести видов 
неоднородности и однородного полупространства неоднородности и однородного полупространства 


На рис. 5 показано изменение вертикального смещения в приповерхностной зоне непре- 
рывно-неоднородного покрытия при воздействии с поверхности постоянного теплового потока в 
пределах круга единичного радиуса, толщина неоднородного покрытия равна радиусу зоны при- 
ложенного воздействия. 





Рис. 5. Смещение поверхности 
непрерывно-неоднородного полупространства 


Величины термомеханических параметров на поверхности покрытия вдвое превышают их 
значения в подложке и линейно убывают по глубине. 
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$.М. АТ2ТКОМТСН, 1.1. КВЕМЕУ 
(Роп Зае Теспткса! Упмег®Ку) 


Тре питепса/ апа!УЙс 5оибоп юЮ те иптихей ахбуттес ргоМетз о! пегтое/азИсйу Гог е сопйпиоис- 
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